Inverse cascade behavior in freely decaying two-dimensional fluid turbulence 
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We present results from an ensemble of 50 runs of two-dimensional hydrodynamic turbulence 
with spatial resolution of 2048^ grid points, and from an ensemble of 10 runs with 4096^ grid points. 
All runs in each ensemble have random initial conditions with same initial integral scale, energy, 
enstrophy, and Reynolds number. When both ensemble- and time-averaged, inverse energy cascade 
behavior is observed, even in the absence of external mechanical forcing: the energy spectrum at 
scales larger than the characteristic scale of the flow follows a k~^^^ law, with negative flux, together 
with a law at smaller scales, and a positive flux of enstrophy. The source of energy for this 
behavior comes from the modal energy around the energy containing scale at t = 0. The results shed 
some light into connections between decaying and forced turbulence, and into recent controversies 
in experimental studies of two-dimensional and magnetohydrodynamic turbulent flows. 

PACS numbers: 47.27.-i,47.32.C-,47.65.-d 
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I. INTRODUCTION 

Forced and freely decaying flows are different, as far 
as long-time energetics is concerned. While forced flows 
can develop both "direct" and "inverse cascades" [ll, [1] 
(transfer of ideally conserved quantities to smaller or 
larger scales with constant flux), freely decaying flows 
can at most display "selective decay" 0] (the faster 
decay of one invariant relative to another, when the sys- 
tem has two or more ideally conserved quantities and one 
of them condenses at large scales in the ideal case) . These 
two processes are in practice related: selective decay re- 
sults from the transfer of one of the invariants to short 
wavelengths where the dissipation coefficients are more 
effective. 

For flows with only direct transfer of ideal invariants, 
e.g., three-dimensional (3D) hydrodynamic flows, the dis- 
tinction between freely decaying and forced flows is not 
that important. It is well known that if time-averaged 
around the time of maximum dissipation, the spectrum 
of the ideal invariant develops in the decaying case the 
same scaling as in forced cases. It is also known that the 
shape of the spectrum during the decay is preserved un- 
der certain conditions Q , resulting in a self-similar decay. 
The concept of "direct cascade" is then often applied to 
these flows indistinctly of the nature of the mechanisms 
sustaining the turbulence. 

However, for flows with inverse transfer of an invari- 
ant, the distinction between freely decaying and forced 
flows remains, and dates back to before the introduction 
of concepts such as selective decay. In two-dimensional 
hydrodynamic flows, on the one hand, Kraichnan consid- 
ered the forced case @ and predicted an inverse cascade 
of energy with a direct cascade of enstrophy. On the 
other hand, Batchelor Q considered the direct transfer 
of enstrophy in decaying flows, obtaining the long dis- 
cussed ^ enstrophy spectrum (see, e.g., [8l-[ll|). The 



reasons for the distinction in this case are clear: as ex- 
plained by Kraichnan, the possibility of an inverse cas- 
cade depends on the relative strengths of the non-linear 
interactions between scales. It is natural then to assume 
that if the input of energy is not sustained for sufficiently 
long times, the inverse cascade cannot develop and that 
in the flow only the direct cascade, with an increase of 
the energy containing scale, will be observed (note that 
the increase of this scale is however distinguishable from 
the three-dimensional hydrodynamic case). 



The main aim of this work is to study under what con- 
ditions, in a decaying flow, behavior associated with an 
inverse cascade (approximately constant negative flux, 
and the corresponding self-similar spectrum with the ap- 
propriate scaling law) arises. The motivation to identify 
such behavior is twofold: On the one hand, phenomena 
observed in solar wind turbulence and space plasmas are 
often loosely associated with an inverse cascade of mag- 
netic helicity [l^ - fl^ . In many of these cases (as, e.g., 
the solar wind) the large scale dynamics can be well ap- 
proximated by a decaying magnetohydrodynamic (MHD) 
flow. Although forced MHD flows develop inverse cas- 
cades (of magnetic helicity in three dimensions [T5| . and 
of the square vector potential in two [11]), based on the 
previous arguments results obtained from simulations or 
theory of forced turbulence do not strictly apply to the 
decaying case. On the other hand, recent experiments of 
decaying soap film flows (see [iTj for a review) showed in 
some cases positive third-order velocity structure func- 
tion |18j . inverse energy transfer for scales larger than 
the injection scale [l^, thickness behaving as a passive 
scalar with Kolmogorov-like scaling [l^] , and in one case 
"^/•^ energy scaling 21|, results that 



even a disputed fc" 
are considered indications of inverse cascades in forced 
two-dimensional flows, but are puzzling in the decaying 
case. 
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II. NUMERICAL PROCEDURE 

The systems discussed above are complex, with small- 
scale plasma effects in the case of space flows [l^] and 
with deviations from two-dimensionality (23j and finite- 
size effects [3,[23l in the case of soap films. These effects 
are important to explain the observations (see, e.g., the 
discussions in [lllifl]). With a more humble objective in 
mind, we take these problems as motivation, and focus 
our attention on differences that may arise in simulations 
and experiments because of the differences in the proce- 
dures followed to obtain turbulent data. To this end, 
we restrict ourselves to a much simpler system: that of 
two-dimensional hydrodynamic flows in the absence of 
external forces, freely decaying in periodic boxes. 

In high resolution numerical simulations, because of 
the computational cost, seldom an ensemble of simula- 
tions with different initial conditions is done. Instead, 
only one simulation is done, and ergodicity is assumed 
in order to obtain statistical quantities by averaging in 
space. Experiments of decaying turbulence are character- 
ized instead by single- or multiple-point measurements of 
quantities downstream for very long times. The multiple- 
point measurements can be at similar stages of the evo- 
lution, or at different stages (e.g., when measuring at 
different distances downstream from the place where tur- 
bulence is generated). The procedure actually results on 
an average over multiple realizations of decaying flows. 

To mimic such a procedure in our simulations, we per- 
formed 50 two-dimensional simulations of the Navier- 
Stokes equations on a regular grid of 2048^ points. Ini- 
tial conditions were a random superposition of harmonic 
modes between wave numbers A: = 18 and 22, with 
the spectrum peaking at ko ~ 20. Viscosity was u = 
2.5 X 10~^ in all these runs, and the box had length 27r. 
The initial r.m.s. velocity Ur.m.s. in all runs was 1, cor- 
responding to a turn-over time of tjsil — Lq /Ur.m.s. ~ 
27r/20 ~ 0.3 [Lq = 27r/A:o)- As a result, the runs only 
differed by their random initial phases. 

To study the effect of scale separation between the ini- 
tial energy-containing scale and the box size, a second set 
of 10 simulations was done with a resolution of 4096^ grid 
points. The initial random excitation with Ur.m.s. ~ 1 
was placed between A: = 30 and 34, with the spectrum 
peaking at fco ~ 32. The initial integral scale of the 
flow was thus decreased by a factor 1.6, and the initial 
turnover time was rjv^ « 0.2. The viscosity in these runs 
was V = 9.9 X 10^^, and as in the previous dataset, the 
runs in these dataset only differed in the choice of initial 
phases. 

The equations were integrated up to t = 6 using a 
parallelized pseudo-spectral code, with the 2/3- rule for 
dealiasing, and second-order Runge-Kutta to evolve in 
time [1^. Neither friction nor forcing were employed. 
Note that although the spatial resolutions considered in 
this study are moderate, the computational cost of com- 
puting the whole dataset (the 50 2048^ runs and the 10 
4096'^ runs) is equivalent to that of computing a 16384^ 
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FIG. 1: (Color online) Time evolution of the energy spectrum 
(a), energy flux (b), and enstrophy flux (c) in a single 2048'^ 
simulation, from t — 0.5 (black) to t = 6 (light gray or light 
blue). Slopes in the energy spectrum are indicated as refer- 
ences. The curves corresponding to t = 0.5 and to t — 6 are 
indicated in all panels by arrows, and the vertical dashed lines 
indicate the initial energy containing wave number ko- In (b) 
and (c), note the displacement to smaller wave numbers of 
the minimum of energy flux, and to larger wave numbers of 
the maximum of enstrophy flux. The inset in (a) shows the 
time evolution of the energy (solid line) and of the enstro- 
phy (dashed line) in this run, with the color changing with 
time following the colors used for the different curves in the 
spectrum and fluxes. 
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FIG. 2: (Color online) Time evolution of E</E (solid) and 
/E (dashed) (resp., ratio of the energy at large and at 
small scales to the total energy) as a function of time for all 
runs in the set of simulations with 2048^ grid points (a), and 
with 4096^ grid points (b). 
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III. RESULTS 

Since the first set has more realizations (50), result- 
ing in more reliable ensemble averages, we focus first on 
this set and compare at the end with the second datasct 
to identify the effect of increasing resolution. Figure 
[T] shows the time evolution of the energy spectrum, of 
the energy flux, and of the enstrophy flux in one of the 
2048'^ simulations, i.e., for a single realization. At large 
wave numbers the energy spectrum develops a direct en- 
strophy cascade range with an energy spectrum steeper 
than ^ but shallower than ^ k~^. Such spectra 
have been reported before, and for details we refer the 
reader to the discussion in [ll| and references therein. 
We simply note that, given the separation of scales in 
this run kmax/ko ~ 34, with kmax = N/3 the maxi- 
mum wavenumber, one does not expect to be able to 
resolve sufficiently the small-scale range. Recent simu- 
lations of forced two-dimensional turbulence, with more 



FIG. 3: (Color online) (a) Time evolution of the integral scale 
for all the 2048^ runs; note the dispersion around the mean 
evolution. The inset shows the time evolution of the energy 
for all runs, which are very close, (b) Time evolution of the 
standard deviation a in the integral scale of all the 2048^ runs. 
The inset shows the standard deviation in the energy in all 
the 2048^ runs as a function of time. 



scale separation between fcg and kmax, indicate a clear 
~ k~^ energy scaling at small scales [23, as expected 
from Batchelor phenomenology (see also[l|, [l^ ) . Since 
this is not the main purpose of this paper, we just want 
to stress that the resulting inertial range is compatible 
with moderate Reynolds numbers runs done in the past. 

We are interested here instead in the spectrum at wave 
numbers smaller than the initial energy containing wave 
number fco. At those wave numbers, the peak of energy 
moves towards smaller wave numbers as time evolves, 
leaving to its left an energy spectrum again consistent 
with the direct enstrophy range. The fluxes arc in agree- 
ment with this picture: the energy flux has a negative 
peak that moves towards smaller wave numbers with 
time, but has no identifiable range with approximate con- 
stancy. On the other hand, the enstrophy flux develops 
a wide quasi-constant range starting at t « 1, and as 
its amplitude decreases with time, this range can still be 
identified. 

The displacement of the peak of the energy spectrum 
towards smaller wave numbers can be associated with the 
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well-known increase of the flow integral length L as the 
energy decays. However, it does not automatically follow 
that the increase of L is the sign of an inverse cascade. 
Indeed, the integral length scale also grows in a three- 
dimensional flow as is well known, due to the preservation 
of large scale correlations in the von Karman-Howarth 
equation [2^, [2^ with spectra ~ /(;+", with n either equal 
to two or four, and also due to the eddy-noise created 
by the beating of two interacting small scales. As the 
peak of the energy spectrum in three-dimensional turbu- 
lence moves to larger scales, the amplitude of the peak 
decreases, being modulated by the ~ A;+" spectra. On 
the other hand, the growth of the integral scale in an 
inverse cascade is accompanied by another phenomenon, 
that of the growth of energy at large scales. 

Unlike freely decaying three-dimensional turbulence, 
the energy at large scales in the two-dimensional runs 
increases substantially as time evolves. This can be seen 
in Fig. [21 that shows the energy at large scales (wavenum- 
bers smaller than fco) 

E<=Y. Eik), (1) 

k=l 

and the energy at small scales (wavenumbers larger than 
fco) 

E>=^E{k), (2) 

k—ko 

normalized by the total energy as a function of time. In 
both runs, the small-scale energy decreases rapidly, while 
the large-scale energy grows. When not normalized by 
the total energy, E'^ still grows and then slowly decays 
as a result of viscous forces (with a slower decay in the 
4096^ run). In the absence of transfer of energy towards 
large scales (i.e., without negative energy flux), £'< would 
only decrease in time. 

The time evolution of the integral scale for all 2048^ 
runs is shown in Fig. [3l In spite of the fact that the runs 
only differ by their initial phases, there is a large dis- 
persion in the time evolution of L, more so than for the 
energy (see the inset of Fig. [3lja)). The actual dispersion 
in the evolution of L and E in the different runs is also 
quantified in Fig. ^h) , by means of the standard devia- 
tion a. For both quantities, a first grows exponentially 
(at very early times) , and later seems to follow an approx- 
imate linear growth with time. From pioneering works on 
predictability of two-dimensional turbulence [30l - l32| (see 
also 33,, ,34J for recent studies) , we can expect differences 
in the initial conditions to grow first at the initial en- 
ergy containing scale (resulting in the early exponential 
phase) , and later to propagate towards larger scales if an 
inverse cascade develops (see [S^l for a numerical study) . 
In that regime, the time it takes for the differences to 
propagate is that of the turnover time, which if a Kol- 
mogorov spectrum is assumed, results in linear growth of 
the error with time [lOjIl^l- Note that the integral scale L 
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FIG. 4; (Color online) Energy spectra for 10 2048^ runs at 
t = 1.5 (no averaging performed here). Slopes are indicated 
as references. The vertical dashed line corresponds to the ini- 
tial energy containing wave number fco. The inset shows the 
same ten spectra at t = 6. Note that the large-scale spectra 
look like ~ k, which could be interpreted as eddy-noise in 
two-dimensions (see [sHi). However, unlike three-dimensional 
turbulence, the amplitude of the peak of the energy spectrum 
is larger than its initial value (at f = 0) at the same wavenum- 
ber. 



is obtained from the energy spectrum weighting the most 
energetic wavenumbers, and as such more deviation can 
be expected than in the case of the total energy. 

Although the standard deviation observed in Fig. [3] 
can be expected in an ensemble, numerical simulations 
often deal with only one of the realizations. What are 
the implications of the dispersion in the integral length 
(and energy) in the different simulations at a given time? 
Figure 21 shows the energy spectrum at t = 1.5 and at 
t = 6 for 10 randomly picked 2048^ runs. At t = 1.5, 
a narrow range seems to emerge at large scales, with a 
slope shallower than ~ k~^, between the energy contain- 
ing scale and the enstrophy range. Exploration of the en- 
ergy flux (not shown) indicates that its ensemble average 
may display a short range with approximate constancy. 
Moreover, at later time (see, e.g., the inset in Fig. H]), 
this shallower range in the energy spectrum increases in 
width and moves to smaller wave numbers. 

Based on this result, we show in Fig. [5] the time- and 
ensemble-averaged energy spectrum, energy flux, and en- 
strophy flux for the set of simulations with 2048^ grid 
points. The average is computed over the 50 simulations 
and from t = 0.5 to < = 6. Now, a wide ^ k~^/'^ spec- 
trum clearly emerges at wave numbers smaller than the 
initial energy-containing wave number, followed by the 
enstrophy range at larger wave numbers. The energy 
flux now also displays a wider range in which it remains 
negative and non-negligible (note that negative energy 
flux corresponds to positive third-order longitudinal ve- 
locity structure function). Moreover, a clear distinction 
can be made between the range with negative energy flux 
and negligible enstrophy flux, and the range with positive 
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FIG. 5: Time- and ensemble-averaged energy spectrum (a), 
energy flux (b), and enstrophy flux (c) over the 50 2048'^ sim- 
ulations and from t — 0.5 to f = 6. Slopes in (a) are indicated 
as references. The vertical dashed lines corresponds to the 
initial energy containing wave number ko- 



enstrophy flux and negligible energy flux. Finally, note 
that, dimensionally, E ~ fco|n|, which is well verified by 
the data; in other words, the observed energy and en- 
strophy fluxes (respectively to large and small scales) are 
comparable, when properly normalized using the charac- 
teristic length scale of the flow at t = 0. 

The same time- and ensemble-averages are shown in 
Fig.Elfor the 10 runs with 4096^ grid points (as the scale 



FIG. 6: Time- and ensemble-averaged energy spectrum (a) 
and energy flux (b) over the 10 4096^ simulations and from 
t — 0.5 to t = 6. Slopes in (a) are indicated as references. 
The vertical dashed lines corresponds to the initial energy 
containing wave number fco- 

separation between initial perturbation and dissipation 
scale is roughly the same in this set of runs as in the 
runs with 2048^ points, the flux of enstrophy is similar 
as the one in Fig. [5] and is not shown). The larger scale 
separation at large scales allows for the development of 
a slightly wider ^ spectrum. 

The results in Figs. [5] and [H] indicate that after timc- 
and ensemble-averaging, well-known features of the in- 
verse energy cascade in forced two-dimensional hydro- 
dynamic turbulence can be identified in decaying flows. 
The ~ energy spectrum emerges as an enve- 

lope to the evolution of the large-scale energy spectrum, 
which instantaneously and for each run only shows a 
clear direct enstrophy range. It is interesting to point 
out that time- and ensemble-averages of freely decay- 
ing two-dimensional turbulence were used before (see 
j36|). There, a spectrum for the enstrophy ~ k (i.e., 
E{k) ~ k~^) was observed at scales larger than the ini- 
tial energy containing scale. Although this indicates a 
growth of energy at large scales, the spectrum is shal- 
lower than the one reported here, probably because in 
that study much later times were considered when study- 
ing the decay (over 100 turnover times). 
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IV. DISCUSSION 

Experiments, as well as observations of space flows, do 
not allow for such a tight control of initial r.m.s. velocities 
and initial injection scales as in the present simulations, 
which kept these quantities the same in all the initial con- 
ditions. It can be conceived that as these initial quanti- 
ties are allowed to vary within a certain percentage, the 
dispersion of the integral scale and of the energy at a 
later time during the decay will become larger. This, in 
a much larger ensemble of runs, may allow for a ^ k~^^^ 
energy spectrum to arise without time averaging, as each 
realization at a given time may correspond to different 
stages of the decay. 

Although in the applications discussed in the introduc- 
tion other effects may be playing a role to explain the ob- 
servations, it is however remarkable that a link between 
the behavior of forced and decaying flows was missed for 
such a long time. We offer a possible explanation for this. 
In a time when computers allow us to study larger and 
larger Reynolds numbers, and many studies are focused 
on achieving the largest possible Reynolds number, the 
present study points at the need to use computer power 
also to perform longer runs, and ensembles of runs (as 
routinely done in climate studies) to be able to compare 
on equal footing with experiments and theory, and to val- 
idate some of the usual assumptions in the field, such as 
ergodicity (in the sense that time- and spatial-averages 
can be seamlessly interchanged when comparing experi- 
ments and simulations), and, as was our focus here, dis- 
tinctions between forced and decaying flows. 

It is worth recalling in this context that the statistical 
mechanics of ideal (unforced) flows, which have a signa- 
ture of condensate of energy at large scales, were the basis 
on which Kraichnan's conjecture of the inverse cascade 
was made, so that it is not necessarily surprising to find 
such a behavior in the decay run as well, since it stems 
from the nature of two-dimensional nonlinear triadic in- 



teractions, and the constraints of conservation. Indeed, 
the triads responsible for the inverse cascade of energy 
are present in freely decaying flows, and whether they 
can act or not is a problem of energetics (i.e., whether 
their amplitude is large enough, and can be sustained 
for sufficiently long times). It is therefore important to 
point out that it has been shown that a ~ k~^/'^ energy 
spectrum is the minimally steep spectrum that is able to 
support an inverse cascade of energy [s^] (i.e., a shallower 
spectrum cannot sustain the inverse cascade). The be- 
havior observed in these freely decaying simulations can 
thus only be observed for a certain period of time, as the 
energy is spread over more wave numbers and the inverse 
energy flux decreases with time (see Fig. [T|) . In the forced 
case, the external forcing supplies the required energy to 
sustain the flux in time and to give a persistent ^ k~^^^ 
energy spectrum. 

Interestingly, recently other problems were found in 
fluid dynamics in which particular triads allowed for tran- 
sient or for sustained inverse cascades. Such is the case, 
e.g., of three-dimensional helical flows restricted to tri- 
adic interactions with only the same sign of helicity in 
all interacting modes [38|. In this later case, the inverse 
cascade was sustained in time instead of being transient, 
as it resulted from an artificial amplification in the sim- 
ulation of particular triads. 

Other similarities between decaying fiows and the 
forced Kraichnan's inverse energy cascade may arise, and 
arc left for future work. In particular, intermittency of 
the large scales, the effect on the structures, and the 
equivalent behavior in MHD can be studied and may be 
of interest for many applications. 
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